An Extended Approximation for the Lowest-lying States in Odd-mass Nuclei 
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An enhanced model, based on the Extended Boson Approximation, for the lowest-lying states in 
odd-mass nuclei is presented. Our approach is built on the Quasiparticle Phonon Model, extending 
it to take into account the ground state correlations due to the action of the Pauli principle more 
accurately than in the conventional theory. The derived interaction strengths between the quasi- 
particles and the phonons in this model depend on the quasiparticle occupation numbers explicitly 
coupling the odd-mass nucleus equations with those of the even-even core. Within this model we 
calculated the transition probabilities in several Te, Xe and Ba isotopes with A«130. 
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I. INTRODUCTION 

Due to its simplicity and numerous successful applica- 
tions the Random Phase Approximation (RPA) [lH3j is 
widely considered as a good first approximation to study 
small fluctuations in atomic nuclei. However, this sim- 
ple model enjoys only a limited success when one needs 
to describe properties of states from the lowest part of 
the spectrum in nuclei remote from the magic configu- 
rations. The Quasi-Boson Approximation (QBA), un- 
derlying the RPA, stimulates a discussion concerning its 
applicability to the problem of correctly taking into ac- 
count the ground state correlations (GSC) in even-even 
nuclei. Numerous improvements of this theory with re- 
spect to adding correlations in the ground states of even- 
even nuclei have been attempted, as for example in 0- 
7]. These enhanced models stem from the disregard of 
the QBA and are related to more precise inclusion of 
the Pauli principle when calculating matrix elements of 
various operators. An enhanced version of this approxi- 
mation, referred to as an Extended RPA (ERPA), which 
was proposed a long time ago [8| and later developed 
in llCf , proved successful in improving the theoretical 
results for most measurable quantities near the nuclear 
ground states as, for example, the transition charge den- 
sities in the interior region. 

In the present work, we follow the ERPA approach, 
extending it to provide a refined version of the Quasipar- 
ticle Phonon Model (QPM) for odd-even nuclei UMM- 
The interaction strengths between the quasiparticles and 
phonons in the presented model depend on the number of 
quasiparticles in the ground state. In this way, the core- 
particle equations couple with the generalized equations 
describing the pairing correlations and the excited vibra- 
tional states of the even-even core, thus forming a large 
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nonlinear system. This model is applicable to open-shell 
spherical and transitional odd-A nuclei where the Pauli 
principle effects are becoming essential as the number of 
nucleons in the unclosed shell increase. 

Our research descends from the studies presented in 
[Til [l5j . There it has been shown that the backward am- 
plitudes in the wave functions of these nuclei play a very 
important role for better agreement with the experimen- 
tally measured spectroscopic factors and the properties 
of the states from the lower part of the energy spectrum. 
The theory in the latter papers is based however on the 
QBA which we intend to improve by taking into account 
the action of the Pauli principle more precisely due to 
the Extended Boson Approximation (EBA) Q. 

Another widely adopted approach to study odd-A 
nuclei presents the Interacting Boson-Fermion Model 
(IBFM).The IBFM, introduced in_El and further ex- 
tended in numerous papers (e.g. [17j|). differs from our 
approach by that the excited states of the even-even core 
nucleus are created by operators of a pure boson nature. 
In the IBFM the core-particle interaction depends on a 
number of free parameters which are usually fitted to 
match the spectrum in the odd-A nucleus. In this re- 
spect the QPM is closer to the interacting shell model 
where this interaction is derived from the dynamics of 
the constituting nucleons. 

This paper first outlines the QPM and its extension 
based on the EBA in even-even nuclei. A comparison 
between the models built on the QBA and EBA is estab- 
lished on the basis of the reduced transition probabilities 
from the ground to the first 2 + state in even-even nu- 
clei with A w 130. In section IIII1 we give the QPM 
theory for odd-even nuclei emphasizing the effect of the 
rcnormalization on the interaction vertices. Calculations 
on the spectroscopic factors and transition probabilities 
between states in some odd-even Te,Xe and Ba nuclei, 
where experimental data are available, are presented in 
section HVl Conclusions are drawn in section fVl 
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FIG. 1. Left panel - quasiparticle occupation numbers pj x 100 in the ground state of 1 Ba within a one-phonon QPM 
and EQPM theory for the sub-shells in the valence shell. Right panel - same as in the left panel but for the quantities 
(jif® PM — n® PM /n® PM ^ x 100, where nj is the number of particles on the level j. The quantities in both panels are plotted 
as a function of the first quadrupole phonon's energy. 



II. EVEN-EVEN NUCLEI 

This section aims to mark the basic building blocks of 
the QPM and its EBA extension (EQPM) for one-phonon 
states. The notations used below are the same as in [10( 
and [3. 

In EQPM one defines the quantities pj , which are pro- 
portional to the quasiparticle occupation numbers in the 
ground state on the level j: 



Pj 



V2j 



where a denotes a quasiparticle ("qp" for short) 



^jm Uj'Qtjffi ( )^ Vjdj_ r 



(i) 



(2) 



The other key constituent of the theory is the phonon 
operators ("ph" for short) denned as 

qU = iE[^' At (^''V)-(-i) A -^-;^oy;A-M)]. 

(3) 

The ground state |) in equation ([1]) is the vacuum state 
for the phonon operators: Qx^il) = 0. 

We study the dynamics of nuclear systems governed by 
the simple Hamiltonian in the form: 



0,p) 



H = E <E(^ - AtKW - : (pt Por ; 

T jm 

-Ij2 kW --( m I m ^)--}- ( 4 ) 



accounting for the nuclear mean field, the pairing 
and the isoscalar multipole-multipole interactions, re- 
spectively. 

If the pairing vibrations are not taken into considera- 
tion then one can obtain jlfj the following modified QPM 
equations describing the states in even-even nuclei: 



1 -T l (2j + l)ll- {1 r 2p ^ Ej - X) )=n (5) 
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(9) 



A / 1 



The emergence of the factors (1 — pjf) takes into ac- 
count the blocking effect due to the Pauli principle and 
requires one to solve the equations above as a system of 
coupled equations. 

The multipole-multipole interaction strengths are 
treated as free parameters in our study. In the numeri- 
cal calculations we kept the quadrupole-quadrupole term 
only because it gives the dominant part of the long-range 
interaction for the determination of the low-lying states' 
properties in the nuclei of interest. One way to fix the 
parameter is to have it to reproduce the energy of 
the first 2 + state (w 2 +) . Since a one-to-one correspon- 
dence between w 2 + and exists, we show most of the 
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FIG. 2. The reduced transition probabilities B(E2\gs — >• 2f) (in units e 2 /m 4 ) in several Te, Xe and Ba isotopes plotted against 
the energy u„+ of the first quadrupole phonon. The solid lines represent the experimental energies and transitions. 



calculated quantities as a function of w 2 + because its val- 
ues are more intuitive and closer to the experiment than 
the corresponding interaction strength values. 

Below we discuss the results obtained within the 
EQPM for the quasiparticle and particle occupation 
numbers as well as for the transition probabilities in even- 
even nuclei. 

The differences between the quasiparticle and parti- 
cle occupation numbers in 130 Ba as a function of the first 
quadrupole phonon's energy within the QPM and EQPM 
are presented in Fig. [T] From this figure we see that the 
smearing of the Fermi surface increases together with the 
strength of the field force (and, correspondingly, ui2 1 de- 
creases). In the right panel we point out that the relative 
difference of the particle occupation numbers calculated 
within the two model variants can reach up to 5 %, as is 
the case for the proton sub-shell 2(1^/2 ■ 

The transition probabilities in odd-even nuclei are di- 
rectly linked to the transition probabilities in their corre- 
sponding even-even cores as it will be discussed in section 
IIVI We therefore perform a comparative study of the re- 
duced transition probabilities B(E2\gs — > 2j) in several 
even-even nuclei within the QPM and EQPM. The tran- 
sition probabilities in the EQPM are given as 



The nuclei presented in Fig. [2] were chosen to be close 
to spherical ones, having E(4f)/ E(2f) < 2.5. From this 
figure we can see that the blocking effect due to the Pauli 
principle exerts a large impact on this measurable quan- 
tity. The obvious superiority of the EQPM in this region 
serves as a motivation to study odd-even systems with a 
core described within the framework of this model. Be- 
sides, the transition charge densities, being related to 
the reduced transition probabilities, were studied in [Io| . 
There it was shown that the application of the EQPM 
leads to a better reproduction of the experimentally mea- 
sured distributions in the nuclear interior. 



III. ODD-EVEN NUCLEI 



B(EX\g.s. -> Aj) = 



jj' 



(10) 



In our treatment the states in odd-even nuclei are 
described as mixed states composed of pure quasipar- 
ticle and quasiparticle x phonon (qpxph) states including 
backward-going amplitudes [T3.ll5| 
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FIG. 3. The matrix elements V(^2d 3 / 2 |^3si/2 x 2+) and W (zy2d 3/ / 2 1 ^2c? 3 / 2 x 2^) in 131 Ba plotted against the energy oj 2 + of 
the first quadrupole phonon in 130 Ba. 



¥ tf (JM) = C Jv a + 



, D oXi{Jv)P] M {JM) - Ej„a JM - ]T F jXi {Jv)P jXi {JM)\), 

jXi jXi 



(ii) 



where P] M (JM) = [atgtj 

jm is the qpxph creation The structure coefficients from (fTTj) and the energies of 



operator and ~ stands for time conjugation, according to the states in the odd-A nucleus can be obtained by mak- 



the convention 5,- 



ing use of the equation of motion method. Conforming 
to the relation (JT]), when calculating the matrix elements, 
we obtain the following generalized eigenvalue problem: 



/ ej V(Jj'X'i') 

V(JjXi) Kj(jXi\j'Xi') W(JjXi) 
W(Jj'X'i') 



\ -W (JjXi) 







-W(Jj'X'i') 
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10 

l-£*(JjXi) 

10 

l-£*{JjXi) 



(C Jv 

—Ej v 
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(12) 



For conciseness, we provide only the leading terms of 
the expressions for the matrix elements. 



V(JjXi) = (\{[a JM) H},P} X i}\) = 
= -^=[l-p j +L*(JjXi)]T(JjXt) ) 



W(JjXi) = (\{[al JM ,H],P] Xi }\) 
\1 - n. 



(13) 



(14) 



Kj(jXi\j'X'i') = -[Ij(jXi\j'X'i') + Ij(j'X'i'\jXi)} = 

(15) 

= 8jj>5 xx >8u>[l - pj + C*(JjXi)](ej + w Xi ) 

- 5 jf 6 xx ,6 u >(l +L{JjXi))\Y,^h)C%{jXi\jXi l ). 



I For the numerical calculations we used a diagonal ap- 

—ejpjipf, [1 - pj + C*(JjXi)]— A{Xi\i)<p\ % > , proximation for C [H]. Below we list the notations en- 
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tering into the matrix elements p^|) - (IT51) : 
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Ij(jXi\fX'i') = (\{PMJM), [H,P+ xli ,(JM)]}\), (16) 
C* JW {jXi\j'X'i') = ttaa; £(1 - p ili O<jV&< 
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X Ai (r)+X Ai ,(r) 



A(Am') 
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E 

r 
r 

E 
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(17) 
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(19) 
(20) 



E 



JJ 



(e 2 



ii' 



O 2 



(21) 



In the limit case pj = 0, the problem in (|12l) is brought 
to the model obtained in [14J . Below we discuss the effect 
of the correlations in the nuclear ground state on the 
behavior of the matrix elements in (|12[) . 

The interaction between the quasiparticles and the 
phonons will naturally become stronger when the smear- 
ing around the Fermi level increases. In Fig|3l the de- 
pendence of sample qp-ph interaction strengths on u> 2 + 
is plotted. The weakening of this interaction within the 
extended model, as compared to the interaction derived 
within the QBA, is getting more salient as the ground 
state correlations increase. It is also worth noting that 
the strengths in the backward direction depend not only 
on the structure of the phonon state |Az) building the 
matrix element W(JjXi) but also on all other phonons 
entering into the sum in the second summand of the rhs of 
equation (1141) . This implies that the higher-lying phonon 
states influence the properties of the states near the 
ground state. We estimated that the contribution of the 
higher-lying phonons to the quantities W(JjXi) can be 
up to 25%. The diagonal matrix elements Kj(jXi\j'X'i') 
exhibit a similar effect due to the second summand of the 
rhs of (|14l) . This sum generates an energy shift, which 
can contribute to the appearance of intruder states in the 
lower part of the energy spectrum, as discussed in detail 

In our previous paper [14j, it was found out that the 
decrease in the energy of the first 2 + state leads to a con- 
siderable growth of the quantities W(JjXi) thus pushing 
the first solution very close to the first qpxph pole. This 
did not allow us to correctly reproduce both the proper- 
ties of the odd-even nucleus and its even-even core using 
the same values for the multipole constants n\ and cor- 
respondingly uj\j . We noticed that the values of u> 2 + in 
the even-even core, which let us reproduce the energies 
of the lowest part of the spectrum in the odd-even nu- 
cleus with reasonable accuracy, were much higher than 
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FIG. 4. Quasiparticle strength distribution (C 2 + E 2 ) of the 
state v2d5/2 in 131 Ba. The quadrupole-quadrupole interac- 
tion strength k< 2 > is kept constant in the calculations within 
the three model versions. 



their experimental counterparts. In this regard the weak- 
ened interaction between the pure qp and qp x ph config- 
urations (Tl~3l) - (TT4T) caused by the quasiparticle blocking 
should yield better agreement between the theory and 
the experiment. 



IV. NUMERICAL RESULTS 

In this section, we present numerical results showing 
the influence of the backward propagating terms on the 
spectroscopic factors and the transition probabilities be- 
tween states in odd-even nuclei using the two approxima- 
tions - QBA and EBA - giving rise to different variants 
of the model. The latter are denoted in a similar way as 
in [l|: 

• QPMJ? - one-phonon model, including Pauli prin- 
ciple corrections (as in [121 ]) 

• QPM_BCK_P - one-phonon model, including back- 
ward amplitudes and Pauli principle corrections (as 
in p3) 

• EQPM_BCK_P - one-phonon model, including 
both backward amplitudes and Pauli principle cor- 
rections having a core described within the EQPM. 

Below we give technical details on the calculations per- 
formed. 

For simplicity, we employed a Wood-Saxon mean field 
with parameters fitted to reproduce the nuclear bind- 
ing energies. In a similar way, the pairing strengths 
G T were obtained to match the odd-even mass differ- 
ences in neighboring nuclei(for details see [3). We in- 
cluded quadrupole phonons only since the quadrupole- 
quadrupole interaction, along with the pairing interac- 
tion, plays a dominant role for the low-lying collective 
states in even-even nuclei, as already pointed out in 
scc|TTl We let the quadrupole strength k^ 2 \ correspond- 
ingly u> 2 + , vary and analyze the dependence of the quan- 
tities of interest on lo 2 + . The phonons' energy cutoff is 
set to 15 MeV. One appealing feature of the QPM and in 
particular of the variant described in this paper is that 
the interaction strengths between the quasiparticles and 
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FIG. 5. Spectroscopic factor for the (d,p) reaction in 131 Ba 
as a function of w„+. 



phonons depend only on the parameters describing their 
internal structure, thereby introducing no extra degrees 
of freedom. 

From the computational perspective solving the alge- 
braic system ©-(11]) is a more challenging task than solv- 
ing the equations of the standard QPM. As an initial 
approximation to the solution of the coupled problem we 
take the solutions obtained from the uncoupled equations 
(i.e. Pj = 0). 

We tested the so-developed approximation on several 
odd-A Te, Xe and Ba isotopes entering into the transi- 
tional region. As it has already been pointed out in secull 
where we investigated the properties of the correspond- 
ing even-even cores, the use of the EQPM improves the 
agreement between the results of the calculations and the 
experimental data significantly. 

First, we head off to investigating the single particle 
components of the wave function. In the model versions 
which take into account the backward amplitudes we 
found a serious depletion of the quasiparticle strengths as 
exemplified in Fig2]for the case of the qp state i/2d^/2- 
We found a similar behavior for the rest of the states 



from the valence shell in all nuclei within the consid- 
ered region. An appropriate experimentally measurable 
quantity to study the single particle strength is the spec- 
troscopic factor (SF) for the (d,p) reaction calculated as: 

S Jv = {CjuUj - Ej v vjf (22) 

From FigEl we see that the value of oj 2 + , at which the 
experimentally measured spectroscopic factor is repro- 
duced, is lower in the case of EQPM_BCK_P than in 
QPM_BCK_P by about 50keV and is therefore closer to 
the energy of the first 2 + state in 130 Ba. In Table HI this 
comparison between the two model versions is extended 
for several nuclei where experimental data are available. 
From there we see a systematic improvement with re- 
spect to w 2 + ranging from 50keV to 150keV in favor of 
EQPM_BCW_P. 

While the spectroscopic factors are influenced mainly 
by the properties of the last, unpaired particle, the elec- 
tric transition probabilities depend strongly on the bulk 
properties of the even-even core. The largest contribu- 
tion to these quantities is due to transitions between pure 
qp and qpxph states represented by the sum in the rhs 
of the following expression: 



TABLE I. Spectroscopic factors for the (d,p) reaction of the 
state 3/2+ in 123 Te, 125 Te, 127 Te and 131 Ba . The second 
column gives the experimental pl| values. The columns 3 
and 4 give the energies u 2 + (in MeV) of the corresponding 
even-even cores calculated within QPM and EQPM at which 
the experimental values of the SF are reproduced. 



Nuclide 


Exp 


«„+, QPM_BCW_P 


w,+ , EQPM_BCW_P 


i23 Te 


0.5 


1.4 


1.3 


125 Te 


0.46 


1.5 


1.3 


i27 Te 


0.38 


1.5 


1.35 


131 Ba 


0.25 


1.05 


1 



B odd (EX; J\V\ -> J2V2) = — \ c JiviCj2V2 e npfj 1 j 2 v Jl j 2 + y^U{J 1 uiJ 2 i'2^i)y/B{E\;g.s. -» XT) ] , (23) 

n -h V / 

where e np is 1 if the unpaired particle is a proton and - if neutron; B(EX;g.s. — > A-j) is the reduced transition 
probability in the corresponding even-even nucleus given by formula (|10[) and 

U{J lVl J2V2\i) = —[Cj 2U2 Dj 2 xi(Jivi) - Ej 2V .,Fj 2Xl {J 1 v 1 )][l + L{J 1 J 2 \i)] + 

(~) Jl - j2+x — [Cj lUl D JlXt (J2V2) - Ej lVl Fj lXi {J 2 v 2 )][l + L{J 2 Jx\i)]. (24) 



In expression (|2"3"|) the terms corresponding to transitions between pure qpxph states have been neglected as being 
small. In systems where the last particle is a neutron we make the approximation: 
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FIG. 6. Same as Fig. [5]but for B(£2|3/2+ 1/2+) in 123 Te, 129 Xe and 131 Ba. 



B dd(EX; J\V\ — > J2V2) — — 5- 



^U{J 1 viJ 2 V2\i)y/B{ L EX;g.a. -> A,) 



-^-U 2 (JiViJ 2 V2Xl)B{E\]g.s. 



Ai), 
(25) 



which stems from the fact that the coefficients 
U{J\V\Jiv-2,\i) are non-negligible for the lowest-lying 
states only and out of these states the transition to the 
first excited state is the strongest. 

The dependence B odd {E2\3/2+ -> 1/2+) = 
B odd (E2\3/2+ -> l/2+)(w 2 +) is plotted in Fig. ©within 
the three model versions. This function shows an almost 
linear behavior in the case of the QPM_P while in the 
calculations which take into account the backward am- 
plitudes a peak emerges. This peak is a result of the 
increased fragmentation in the latter pair of model ver- 
sions (cf. FigHJ which contributes to the enhanced values 
of the coefficients U{J\Vi JiVi\%)- As a result, the maxi- 
mum value of the presented transition probabilities in the 
EQPM_BCK_P and QPM_BCK_P is about three times as 
large as the maximum value obtained within the QPM_P 
bringing us closer to the experimental values. It is also 
worth noting that the values of w 2 + , which correspond to 
the peak values obtained within the EQPM_BCK_P, are 
about lOOkeV lower than in the QPM_BCK_P. We there- 
fore conclude that the effect of the renormalization yields 
better results with respect to the experimentally mea- 
sured energy of the corresponding even-even core though 
it is still rather higher from it. 



V. SUMMARY AND OUTLOOK 

In this work we extended the model presented in [l4| 
by taking into account the blocking effect due to the 
Pauli principle, following the approach prescribed in 8] 
and [10j . Renormalized quasiparticle-phonon interaction 
strengths in both the forward and backward directions 
have been derived. Numerical calculations on the spec- 
troscopic factors and transition probabilities in several 
Te, Xe and Ba isotopes have been performed using a 
Wood-Saxon potential well and residual interaction of a 
pairing+quadrupole type. The results indicate an overall 
improved description of these experimentally measured 
quantities due to the weakened quasiparticle-phonon in- 
teraction strengths. 

However, despite the adoption of this elaborate ap- 
proximation, further improvements of the theory towards 
weakening of the qp-ph interaction could resolve some of 
the existing discrepancies with the experiment. Some 
steps in this direction would be the inclusion of higher 
multipolarities, the use of multi-phonon configurations 
and the development of a more elaborate approach to ac- 
count for further correlation effects. We finally conclude 
that our understanding of the properties of the lowest- 
lying states in relatively stable odd-A nuclei still lacks 
the desired accuracy [19 ] . 
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